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Abstract. We develop a simple tensorial contraction method to obtain analytical formula for X-ray resonant 
magnetic scattering. We apply the method considering first electric dipole-dipole and electric quadrupole- 
quadrupole scattering in the isolated atom approximation and compare the results with previous works. 
Then we apply the method to derive phenomenological original formulas which account also for non- 
spherical systems and for dipole-quadrupole mixing. 

PACS. 61.10.-i X-ray diffraction and scattering 33.55.Fi Other magnetooptical and electrooptical effects 



S 1 INTRODUCTION Simple formula for RMXS amplitude can be found in 

^ ■ the literature with the approximation of spherical atoms [3] • 

o ■ 

-ray magnetic scattering exhibits huge enhancement at 

> '■ 

■ resonances corresponding to localized electronic states |ll2| . 

X 

^ Because of its coherent nature and dependence on polar- ^hey are composed of a sum of terms, each term be- 

isation. X-ray resonant magnetic scattering (XRMS) can jj^g the product of a geometrical expression, containing 

be used to determine both spatial distribution and orien- polarisations, multiplied by a complex valued function of 

tation of matter magnetisation. energy. In a phenomenological analysys, where one is in- 

To design experiments and analyse data, there is a terested in separating the different contributions to the 

strong need for simple analytical formula which directly spectra, one can take theses complex valued function as 

relate the experimental geometry and sample orientation free parameters, under the constraint that dispersion re- 

to the measured scattered intensity. lation are satisfied. 
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We derive a simple and understandable theoretical fra- above equation by the following definition where one has 
mework within which we recover previous XRMS formula factored out a icr^X/V factor ( rg is the classical radius 
and then we extend such formula to non-spherical systems, of the electron ) : 

tek^e'k' = ic/VroXF^k^^>k' (4) 



2 introduction to scattering factors 



The interaction between matter and a photon described 



or 



-, ( 0|e-^'''-''p • e'|7i ) (n|e*-''p • e|0) 
by a wave vector k and polarisation e, is written, dis- p^^ _ _\ /__ 

~^ fim ^-^ a; — a;,, + 

n 

carding elastic Thomson scattering ( the A {r) term in 

the Schroedinger equation which becomes dominant off- ^^iis formula is the starting point of our analysys. To 
resonance), and discarding also the spin-magnetic field in- analyse such formula we expand exp{ik ■ r)p ■ e to first 
teraction fH , as : order, and rewrite the electron momentum p as ^[H, r] : 

H,nt = (1) exp{ik ■r)p-e~{l + ik-ry-^[H,r-t] (6) 



{a\^,exp{-i\<. ■ r) + afe,,ea;p(ik • r))— p • e =[H,r-e]+ i[H, k ■ rr ■ e]/2 



ujV J ' ' cm 

(2) where the exponential has been expanded up to first order 

where r is the electron coordinate, p the kinetic mo- and H is the system Hamiltonian. Such substitution leads 

ment, a^k,e, o,k,e are creation annihilation photon opera- to : 

tors, Lo = kc, and V is the space volume. j-, 

The photon-photon resonant scattering is a second or- ^ ^2 ^o|(e' r-^fc' re r)|n) (n\(e-r+^k-re-r)\o) 

h uj—uJn+iS 

der perturbative process whose resonating amplitude per 

(7) 

unit of time, for initial e. A; and final e',/c' states, is ob- 

tained from matrix elements as : we are going to describe atoms as spherical entities 

/„, ,1,' r n \ / I y plus some kind of distortion it is interesting to consider 



hujm?V ^-^ bj — ijjn + id the spherical group and how the angular part of the dif- 

(3) 

ferent terms of Eq. |^ decompose on the irreducible rep- 

where < n| is a complete set of eigenstates of matter, , , . r 1 mi 1 1 , 

resentations of such a group, ihe scalar product e • r can 

< 01 being the initial state, hujn is the energy difference , „ 1 , r 1 , 1 • 1 i 

be rewritten as a sum of products of rank 1 spherical har- 

between < n\ and < 01 and we have retained only the 

monies tensors components : 

resonating denominator. In the literature one always finds g=+i 

e . r = ^ el*rl (8) 
the scattering amplitude Fcfc^e'fc' which is related to the o=-i 
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where, given a cartesian vector A, its spherical tensor com- while in section 21 we consider a non-spherical perturba- 
ponents are tion term and derive a phenomenological expression for 

the magnetic scattering amplitude. 

Aq = Az 
Al = -{A, + iAy)/^i2) 
AU^ {A^-iAy)/^{2). 

(^gj 3 derivation of formula for the spherical case 
with a magnetic field perturbation 



The product k • re • r is the sum of a constant contribu- 
tion ( rank 0) and a sum of products of rank 2 spherical 
harmonics tensors: 

q=+2 



In the case of a spherical atom perturbed by a magnetic 
(k • r)(e • r) = o ' ^)('" ' ^' 5Z ^^^)q^'''^)<i exchange field the final scattering amplitude expression 



q=-2 

where, given two cartesian vectors A and B , the sym- 
bols (AB)q denote the rank 2 spherical symmetric tensor 
components of their product which are : 



does not mix dipolar with quadrupolar terms because of 
parity conservation. We take the angular moment quan- 
tisation axis 1^, along the magnetic field. Terms having 
different quantisation number q will not mix. The expres- 

i2A B — A B — A B ) 

y/Q ^ ^ ^ ^ y y' sion for scattering then assumes the following simple form 

1 

2* 

1 

^ q=l q=2 

(11) F.k^e'k' = J2 P^^i'l*4 + E F2,qike);ike), (12) 

q=-l q=-2 



(AB)±i = T:^iA,B, + A,B,, ± i{AyB, + A,By)) 



(AB)±2 = -{A^B^ - AyBy ± i{A^By + AyB^)) 



We see from Eq. ^| that spherical tensorial compo- 
nents of the tensor product k e are coupled to same rank 
components of the product r (g) r to form a scalar. As the 
product r §5 r has no antisymmetric component the rank 1 
components of the product k (g) e do not appear in Eq. IIUI 

The transformations expressed by Eas. I9I1I1 are used 
to simplify the scattering factor expression by choosing 
the angular moment quantisation axis in the most ap- 
propriate way. In the next section we consider a spheri- ^2,9 = ^2 + 9-^2 + Q^^i + Q'^^i + 9^^^ 
cal system pertubed by a magnetic exchange interaction (13) 



In the spherically symmetric case (no magnetic field) 
the Fi q and F2.q are independent of q. The introduction 
of a magnetic field introduces a q dependence in the scat- 
tering factors that can be expanded as a polynomial in 
the quantisation number q : 

Fi^y = + qFl + (z^f 2 
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Where the coefficients are given by : For dipolar scattering we have therefore 

= (2i^i,o - - i^i,-i)/2 



such formula can be rearranged in the Hill & Mc Morrow 
[3| form : 



-^2" — ^2,0 

F2 = (^2,-2 ^ -^2,2 + 8i^2,l ~ 8-F2,-l)/12 

Fl' = (I6F24 + 16F2,_i - F2._2 - ^^2,2 - 30F2,o)/24 

Fl' = (^^2,2 - ^^2,-2 + 2^^2,-1 - 2F24)/12 

F|' = (6^^2,0 - F2,2 - 2^2,-2 - 4^^2,1 - 4F2,_i)/24 Ffj^^l = e' ■ e« - F[^) ^i■ (e' x e)F[^ + {e' ■ i){e ■ 



(14) 



(17) 



In the non magnetic case all the terms are zero except For the quadrupolar term the same procedures apply 

the zero order ones. with some more bookkeeping for the various terms that 

We can substitute q in Eq. [T^lby |.L where L is the arise when we apply the i^x operator on the right terms, 

angular moment operator. Going back to cartesian space, In details : 
Eq. Elis then written as 

F,k^,,k, - C ( e' ^ {tix^F^^e] + 



(jCx)°e(g) k = e(g) k 

C [ k'(g)e'^ (i|x)"i^f k® e ) /2 {i£_x)h<^k = i{i x e) k + ® (^~ x k) 



\ ri=0 
n=4 



\ ri=0 / 

(15) -{£_xfe(^k = (g) k- 2(^ X e) (g) (^ X k) + e(g) ki 
iifx (g) k = X e) (g) k + 3ic± {g (£ x k) 

where the symbol C means the sum of all possible contrac- 
tions of the two vectors e' and k' with the expression to + 3i(C x e) « k_L + j(^ x k) g) e 
their right. A contraction is realised by coupling by pairs (Cx )■*£ «) k = «) k - 8(C X e) ® (^ X k) 
the 2A^ vectors which enter the expression , where 2N is _l_ g^^ (g) k^ + k^ (g e 

the sum of the tensors ranks. Each vectors couple is then /, „s 

(18) 

contracted, the result being the scalar product, and the 
N scalar products are multiplied together to give the final 

result. The contraction of a tensor of defined rank with Each of the terms of the above equation must be multi- 

itself must not be considered. plied by its own and contracted with e'k'. The result 
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is : The tensor T must be invariant under all operations 



of the point symmetry group of the system|2]. 

We treat perturbatively the crystal-field correction to 
scattering amplitude considering a process, in one-particle 
approximation, where an electron is promoted from a closed 
shell state, denoted by \nglg^; Jajz), to the intermedi- 
+3iF2^S{e ■ z e- z z-{k' X k))/2 ^^^^^^^ InJamaz) \ ^cr) and Inbhmbz) \ ^cr) of two opens 

+6F2'^e' ■ z e - zk' ■ zk- z shells (ria, L) and {rib, h)- The perturbation T has matrix 

elements which mix the two shells. 



pguadrupolar ^ ^ 2F^2 ^ gF^ )5(e' • £ k' • k)/4 

- 4Ff • (e' x e) k' • k)/2 
+ _ 7F^'')S{e' -ze-zk' ■ k)/2 

+ (-2F^2 _ SF^'^)S{z • (e' X e) z • (k' x k))/4 



(19) 

where the expression S{t) means the symmetrised ex- 
pression formed by the term t plus the other three derived 



We write the initial one-electron state, discarding some 
unessential labels, as: 



terms that one obtains swapping e' with k' and/or e with i 7 ■ \ _ \^ f ■ M \ 17 ■ \ (ot^ 

I'JaJz) — / ^ C[Jz, (T )\(T > \lg, Jz — (T > (^/ij 

k. Our formula for quadrupolar scattering is in agreement 

with HiU & McMorrow onejSl. '^^"''^ ""^^^'^^ Clebsh-Gordon coefficient 

C{lg,jz - <jA/2,(J;Ja,jz)- 



4 Non spherical case. 

In this section we develop a phenomenological formal- 
ism for resonant magnetic scattering in the non-spherical 
case. Our starting point is crystal-field theory where the 
non sphericity of the atomic environment is represented 
by a one-particle mean field potential which is added to 

the atomic Hamiltonian. Such an approach, pioneered by '^frng^+a- a) < ^51 "^gz 1^6* \h, iTigz + q' > x 

rrig- <7 qq' 

the works of Bequerel|5|, BethejB], Kramers|3 and van 



The electron-photon interaction is represented by the 
tensors P' ~ P{e' ,k') for the outgoing photon, and P for 
the incoming one. 

The scattering amplitude at first order in T takes the 
form : 



VleckjH], was applied to the calculation of x-ray absorp- 
tion spectra and scattering factors for the first time by 
van der Laan and Thole Pj. 

The crystal field T is given by a superposition of spher- 
ical tensors : 



fb{mg,+q'.a) < h, rUgz + q'\T\la, mgz + q> X 
fa(mg^+q,a) < ia,1Tlgz + q\Pa\lg , "m-gz > 

+ [a ^ 6] (22) 

In this expression the factor /^(to, a), or faifn, a), con- 
tains implicitely the electron propagator for the spherical 



T = Y.^i,qT\ (20) 

; g atom and accounts also for the orbital occupancies. 
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c^irriz + a, a) for L3 


rriz = 1 rriz = = —1 


a = 1/2 


1 2/3 1/3 


a = -1/2 


1/3 2/3 1 



Table 1. values of c{jz,a) at the La resonance 

Starting from this expression we can arrive at a con- 
tracted form which contains, beside polarisation vectors 
and the magnetisation axis, the T tensor representing the 
crystal field, and the spherical tensors Tm^ representing 
the core state. The perturbative process concerns inter- 
mediate levels of well defined angular moment la and li,. 
This imposes a restrictions on the possible contractions : 
the transition due to a P tensor from the L = Ig ground 
state to the L — In excited levels implies that exactly 
{Ig +rank(P) — /„)/2 contractions must be taken between 
the tensor and the inital ground state. Our one-particle 
approximation neglects the energy spread of the interme- 
diate states due to many-body effect and is therefore sim- 
ilar to the fast-collision approximation employed by Marri 
and Carra for the case of dipole-quadrupole scattering in 
a magnetoelectric crystals 

We develop /a(<Z, cr) and fbiq, cr) in powers of L^, as in 
the previous section : 



where the expression is symmetrised by the following 
substitutions: 

symm. = [a b] 

and the symbol |^ (|^) following a polarisation tensor 
means that the number of contractions of such tensor 
with the preceding (following) ground level tensorial ob- 
ject is constrained as discussed above. The (i^mx) opera- 
tor operates on all the object at its right ( remember that 
x)AB = x)A) B + A {{iU X . 

A formal derivation of these contraction rules is given 
in appendix A. 

As an example we give a specialised expression for the 
is edge. We give in tabled the coefficients c^(j2,cr). It is 
useful to notice that for a given — jz — <y the sum of two 
coefficients of opposite spin is constant and equal to 4/3, 
while the difference goes from —2/3 to 2/3 with a linear 
dependence on rriz . This behaviour does not depend on the 
particular edge that we have considered. In general, for an 
edge e, we can always write : {mz +(j,a) = l^, + 2se airiz 
where /g and Sg are constants which depend on the edge. 
For example, at L^, we have If. = 2/3 and Sg = 1/3, while 
at L2 le — 1/3 and Sg — —1/6. 

It is then useful to rewrite equation 1241 as : 



=^/„,,<z"->^/„,,L^ (23) 

"1/2 

The scattering amplitude is then written as : e-ifrpig* p/*\^ / t ^\ni>T''^ (if: vV»" P I^t' 



^ 3(-^n6,l/2/n„,l/2 + /n6,-l/2/na,-l/'^' ^ 



1/2) ^ 



fnb,a-fna:<j'^{mg^+a,cr) +'^e 3 (-^rib , ,1/2 ~ , - l/2^na , - 1/ 

+symm. (24) +symm. (25) 
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In the second term we have converted the quantization The second class is given by the remaining terms which 
number into the operator (i^,„x) at the left of the core- contain the factor v- ((j(^,„x)h), where v is a vector which 
hole orbitals. Once contracted, the above formula gives in enters the composition of an intermediate tensors and h 
the general case a complicated expression. enters the composition of the core-hole tensor. This factor 

However for terms which are zero and first order in can be rewritten as -((i^„iX)v) • h and once again the 
(iS.m'x), the sum over core- hole orbitals can be simplified sum over the core-hole orbitals disappear from the final 
in an elegant way. expression. 

The zero order terms describe anysotropy of the crystal The general expression for scattering amplitude, up to 
field. The first order one describe crystal field induced first order in (i^mx) is therefore : 
corrections to the magnetic scattering. 

Concerning the zero order terms, we notice that such na+nb=i 
sum looks like a trace. In other words if, in a given contrac- n„,ni,=o 

tion diagram, a vector A is contracted with the A;et ground ^i-^b (*^mX) Tm^{iS,m^) Pa) 

orbital, and another vector B is contracted with the same +symm. (26) 

bra core orbital the result is equivalent to a contraction of 

where the an^,nb depends linearly on the / electron 

A with B because the sum runs over a complete basis of 

propagators and the exact linearity coefficients can be 

an irreducible supbspace. 

found working out the contraction diagrams. 

Concerning the first order terms, («^mX) operates ei- 

For higher order terms this simplification is no more 

ther on the corc-holc orbital, or on the intermediate ten- 

appliablc in the same easy way. One could still write con- 

sors ( P, P' , and crystal field tensor T). 

tractions where the core-hole tensor disappear but the 

The terms where {i^m x ) operates on the intermedi- 

formula would be complicate by the fact that when one 

ate tensors can be simplified as in the zero order case 

rewrites v- ((i^„iX)^h) as ((jCmX)^v) • h, the (i^mX)^ op- 
described above. 

erator cannot be factored out of the tensor composed by 

Now we consider the terms where one {i^m x ) operates 

V because 

on the core-hole orbital. The contraction diagrams arising 

from these terms can be split in two classes. («CmX)^(vi V2) 7^ ((iCmX)^vi) V2 -I- Vi ((i^mX)^V2) 

In the first class (i^mx) operates on a core-hole vector (27) 
and the result is contracted with another core-hole vector. The core orbital disappears from the final expression 
This class gives zero contribution because (i^mX) is an because, in our approximation, our L3 core state in Eq.[^ 
antisymmetric operator. is not coupled to the valence orbitals. Such an approxima- 
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tion is valid for the hard x-rays domain where the core-hole The quadrupole-quadrupolc scattering, discarding those 
is deep. terms that can be merged with the spherical formulas, has 

For the sake of simplicity, in the following, examples the form : 
we limit ourself to consider only zero and first order con- 



tributions. 



{{{{ 

(k'•z)(k•z)e'•(^^„xe) 

4.1 Application :Magnetic diffraction amplitude for _l_ [j| x e k] } + [k <-> e] } + [[k' e']] } 

spiral antiferromagnetic holmium / /i ^ 

- [k, e ^ k , e ] } 

We consider the case of holmium. The holmium crystal (30) 
has a hep structure where the atoms are embedded in a 

local D-^h symmetry environment ^^^^ term also contains two scalar products with the z 

Taking a cartesian x, y, z frame with x along the hep axis adding complexity to the amplitude dependence on 

Q g^jj-jg . experimental geometry. 

2 / 2 , 2 , 2nn I , / 3 q 2n , /ooN Thesc tcrms contribute, in our system, to the ampli- 
T = t2{6z - [X +y + z )) ± tslx"^ - 3xy ) + (28) 

tude of the 2n ± q Bragg orders diffraction peaks, n and 

where the first omitted term of the serie is a rank 4 com- 

q being the antiferromagnetic wavevectors. In an experi- 

ponent. The ± signs alternate from one ab plane to the 

ment which measures several 2n±q Bragg order diffraction 

other. 

peaks, the incidence angle may vary considerably from one 

The contribution of the (2z^ ~ x^ — y^) , is centrosym- 

order to the others, and therefore we predict that the in- 
metric. So we consider separately the dipole-dipole and 

elusion of our correction could improve the fit. 

the quadrupole-quadrupolc scattering. Within the frame- 

The scattering factor contribution from the alternating 

work of our simplifying assumption f equation I2()|) wc find 

term ±(x^ — 3xy^) gives diffraction peaks at 2n + 1 + mq 

that the dipole-dipole scattering correction is proportional 

orders, n and m being integers and q being the antifer- 

to 

romagnetic wavevector. To simplify the following treate- 

3e' • z {i(m X e) • z - e' • {i(m x e) - [e' ^ e] (29) , . , , . 

ment, and to give an example oi a term contributing to 

The second term of this expression merges with the previ- the 2n + l + q scattering, we consider magnetisation in the 

ously found form for the scattering in the spherical approx- xy plane and we limit here to the case where, in equation 

imation. But the first term contains two scalar products 1261 only the ria = 1, fib = term is considered. We ap- 

with the z axis and this adds complexity to the amplitude ply the contraction rules to formula 1261 where P' and P 

dependence on experimental geometry. are e and k (g) e tensors respectively. In this we get case a 
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scattering amplitude proportional to : 

(( 

— 6i(e' ■ x){e ■ ■ (x x k) 

+6i{e' • x)(e • y)^™ • (y x /c) 
+6i{e' • j/)(e • x)i.m ■ {y x k) 
+6i{e' ■y){e- y)im ■ {x x k) 
) + [e^k] 
) + [/fc ^ fc', e ^ e'] 



(31) 



occurring at 2n + 1 + mq Bragg orders. Our work com- 
pletes the Hill & McMorrow one|3]. 
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5 Conclusions 

We have established a contraction method by which we ob- 
tain phenomenological analytical expressions for the scat- 
tering factors in terms of scalar and vectors products of 
the polarisation vectors, of the magnetisation axis and of 
the vectors defining the crystal field tensor. Using a per- 
turbative approach we have been able to go beyond the 
usual 502 approximation ^ considering simultaneaously 
the magnetisation and a general crystal field tensor. Our 
method is a viable diagrammatic technique based on con- 
traction between vectors. With our method useful formula 
can be straightforwardly obtained in terms of the polari- 
sation vectors, the crystal field vectors and the magneti- 
sation direction. 

We have applied the method to the case of holmium, 
giving the crystal field correction for the quadrupole-quadru- 
pole scattering and predicting dipole-quadrupole peaks 
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7 Appendix A: derivation of contraction rules 



We show in this section how an angular integration of a 
polynomial over a unit sphere surface can be expressed as 
a sum of contraction diagrams. 

The polinomial integration that we consider is : 



» 27V 

" s 



(32) 



where 2N is the polynomial order, are vectors, and bfn 
is a unit vector. By multiplying and diving this expression 
by J exp{—r'^)r'^^'^'^dr we transform it into an integral 
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over the whole space: 5. J. Becquerel: Z. Physik, 58: 205 (1929) 

2N 



J s 



6. H. Bethe: Ann Physik, [5], 3: 133 (1929) 

7. H. A. Kramers: Proc. Acad. Sci. Amsterdam, 33: 953 (1930) 
y"(]^ a, • Y)exp{~r'^)dv^ /{r{N + 3/2)/2) = g. J. H. Van Vleck: "Theory of Magnetic and Electric Suscep- 



iS 

2N 



2N I, tibilities," Oxford University Press, Oxford and New York, 

1932. 



2N 

\{&^ ■ 5v( / exp{-r'^ + V • v)dv^)\^=n/ {r{N + 3/2)/2) 



2N 

n„ / 2 I A\\ Ti o in\ I -n { T\T I o /o^ 9. Ivan Marri and Paolo Carra : Phys. Rev. B 69, 113101 

• dv ea:p(v /4)|„=o^(3/2)/r(iV+3/2) ' ^ 

(2004) 

(33) 

In this expression, when w = only the terms fully con- 
tracted remain because if a/ • takes down a v/2 factor 
from the exponent, another a,„ • must be used to derive 
this v/2 factor which otherwise gives zero contribution, 
and this gives the scalar product a; • a^. 

The intermediate sum which appear in equation I22L 
can be considered as a projector pi^ over a tensorial space 
of definite rank la- The projection over definite rank space 
is obtained through the |^ and |^ constraint. When ap- 
plying contraction rules to expressions formed by tensor 
of a defined rank, one must bear in mind that, by defi- 
nition, a defined rank tensor gives zero when contracted 
with itself. 
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